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Effective field tlieories that describes the dynamics of a conserved U(l) current in terms of "hy- 
drodynamic" degrees of freedom of topological phases in condensed matter are discussed in general 
dimension _D = d+1 using the functional bosonization technique. For non-interacting topological in- 
sulators (superconductors) with a conserved U(l) charge and characterized by an integer topological 
invariant [more specifically, they are topological insulators in the complex symmetry classes (class A 
and AIII) and in the "primary series" of topological insulators in the eight real symmetry classes], 
we derive the BF-type topological field theories supplemented with the Chern-Simons (when D is 
odd) or the 9-tevm (when D is even). For topological insulators characterized by a Z2 topological 
invariant (the first and second descendants of the primary series), their topological field theories 
are obtained by dimensional reduction. Building on this effective field theory description for non- 
interacting topological phases, we also discuss, following the spirit of the parton construction of the 
fractional quantum Hall effect by Block and Wen, the putative "fractional" topological insulators 
and their possible effective field theories, and use them to determine the physical properties of these 
non-trivial quantum phases. 
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Topological phases are fully quantum mechanical 
states of matter which are not characterized by spon- 
taneous breaking of a global symmetry of the quan- 
tum mechanical system. While gapped in the bulk, 
quite often these phases have gapless excitations at their 
boundary, signaling the highly entangled nature of their 
ground states. Since the discovery of the integer quan- 
tum Hall effect|-|(IQHE) and the factional quantum Hall 
effect (FQHE),EJ the list of topological phases in nature 
has been expanded, in particular, by the recent discov- 
ery of time-reversal symmetric topological insulators in 
two and three dimensions |(2£1| and 3D) in systems with 
strong spin-orbit coupling pilii and the identification of 
■^Hc-B phasj&^s a topological superconductor (actually a 
superfluid) .E3 Unlike the IQHE, the topological character 
of these topological insulators and superconductors (i.e., 
the stable gapless edge or surface modes) is protected by 
time- reversal symmetry (TRS). The presence or absence 
of a topological distinction among gapped phases for a 
given set of symmetries and for given spatial dimensions 
can be studied systematically, and is summarized in the 
"periodic table" of topological insulatpps. and supercon- 
ductors for non- interacting fermions.lijt3 

One of the defining properties of topological insulators 
is their characteristic response to external electromag- 
netic fields. A classic example for this is the quantum 
Hall effect (QHE) which is characterized by the non- 
zero, quantized value of the Hall conductance. The three- 
dimensional time-reversal symmetric topological insula- 
tor is characterized by non-zero magnetoelectric polar- 
izablility, which, in the presence ofLlime-reversal sym- 
metry, also takes a quantized valucliiJ Correspondingly, 
these responses are described by a field theory supple- 
mented with a term of topological origin, such as the 
Chern-Simons (CS) term, or the axion term. For a wide 
class of topological insulators and superconductors, such 
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topological resmniR-theories have been studied and care- 
fully classified. M'tjLJ When it exists, such topologically 
protected response is a powerful way to characterize a 
topological phase, as it is not limited to non-interacting 
topological phases. 

There is, however, a subtle but profound distinction 
between a topological phase and a phase with a topolog- 
ically protected response. For instance, the IQHE is a 
phase of a (generally interacting) two-dimensional elec- 
tron gas in a strong magnetic field whose Hall conduc- 
tance is precisely defined, and protected, by a topolog- 
ical invariant .E3 Similarly its edge states are also pro- 
tected, in this case by their chiral nature. However, the 
FQH states are topological fluids. In addition to having 
a topologically protected fractional quantum Hall con- 
ductance, these fluids are characterized by having non- 
trivial ground state degeneracies on topologically non- 
trivial surfaces, excitations that carry fractional charge 
and fractional statistics (Abelian or non-Abelian), and 
by a set of edge states with generally complex proper- 
ties. The universal properties of such topological phases 
are given in terms of effective low-enerm hvdrodynamic 
topological field theories of these fluids .Ejl3 

While powerful, the effective flcld theories of the re- 
sponse of a system to an external local probe should be 
distinguished from internal or hydrodynamic topological 
fleld theories describing the global properties and the ex- 
citations of lopological phases, as in the example of the 
FQH fluids.Ej Just as fluid dynamics is an efficient de- 
scription of a collection of a macroscopic number of par- 
ticles, the dynamics of topological insulators may be well 
described in terms of collective degrees of freedom, rather 
than relying on microscopic electrons (quasi-particles) . 
In fact, such a picture was developed for the QHE, and we 
have a good understanding of a QH system as a droplet of 
electron liquid. The hydrodynamic topological field the- 
ory description of the quantum Hall droplet is given by 
the Chern-Simons (CS) gauge theory. Once such effec- 
tive description of the low-energy physics is established, 
it is likely to be robust against interactions, and has a 
wider range of applicability than the non-interacting mi- 
croscopic system. The purpose of this paper is to develop 
such a hydrodynamic effective field theory that is capable 
of describing collective excitations in topological insula- 
tors in general dimension. For previous studies, see, for 
example, Rcfs. |4|-|^ 

The CS description of the two-dimensional QHE can 
be derived, for example, from composite particle theories. 
In the composite boson thcory,c3 one first attaches a unit 
fiux to each electron to make it a boson. For the com- 
pletely filled Landau Level, on average this will cancel the 
external magnetic field, and thus we have a composite bo- 
son system which is at zero magnetic field, but interacting 
with a statistical CS gauge field. A subsequent duality 
transformation delivers the CS theory which gives us a 
hydrodynamical description of the filled Landau le 
Similarly, by using the composite fermion picture^ 
the same CS theory can also be derived.EJ 



The flux attachment concept is very successful to de- 
rive the hydrodynamic theories of the Laughlin and Jain 
states of the FQHE. However, it only works in two di- 
mensions since it is based on the existence of the Chern- 
Simons gauge theory and a connection with the braiding 
group. In 2d the excitations (or vortices) can carry the 
quantum numbers of representations of the braid group 
which allows for fractional statistics. However, the con- 
cept of braiding is only meaningful in two spatial dimen- 
sions since the braiding of particle worldlines is not topo- 
logically stable in other spatial dimensions. In addition, 
for technical reasons, flux attachment is somewhat cum- 
bersome to apply even for to 2d lattice systems (see e.g. 
Rcf. Thus, an alternative procedure to derive effec- 
tive hydrodynamic fleld theories of topological phases is 
desirable. 

One alternative approach for deriving effective field 
theories is bosonization. This method is closely related 
to the well-known boson- fermion correspondence in 1 4-1 
dimensions which is an exact operator correspondence 
between a theory of massless (Dirac) fermions and a the- 
ory of a massless relativistic scalar field. We will show 
here that one can use the approach known as "functional 
bosonization" to derive effective hydrodynamic field the- 
ories for general topological phases both in two and three 
dimensions. Although for spatial dimensions d > 1 the 
functional bosonization approach is not an exact map- 
ping, it is nevertheless useful to derive effective low en- 
ergy theories in massive phases, which is what we are 
after in this work. In particular, in this paper, we will 
apply functional bosonization to the (3+l)d topological 
insulator to derive the an effective topological hydrody- 
namic theory of the form of the BE theory in 3-1-1 space- 
time dimensions. 

This paper is organized as follows: In Sec. we start 
by applying functional bosonization to the topological 
insulators in the prirnary aorint; of the periodic table of 
topological insulators £jtj We derive BF-type topologi- 
cal field theories supplemented with a term of topologi- 
cal origin. Next, the physical picture suggested by func- 
tional bosonization is described in Sec. HI, making some 



comparison with flux attachment, and composite parti- 
cle theories. In Sec. [V, we discuss hydrodynamic topo- 



logical field theory descriptions of Z2 topological insula- 
tors, including the three-dimensional time-reversal sym- 
metric topological insulators, and the two-dimensional 
quantum spin Hall effect. They are obtained from the 
BF-type topological field theories of the primary series 
by dimensional reduction. N.B.: We will use D to rep- 
resent space-time dimensions and d to represent space 
dimensions. Thus D = 2 is 1 + 1 dimensional space-time. 

In Sec. we explore possibility of topological phases 
that arise because of strong interactions, e.g., "frac- 
tional" topological insulators. Starting from the BE 
theories for non-interacting topological insulators, and, 
foUoiwM the parton approach pioneered by Blok and 
WcnEj'til in the FQHE by implementing strong corre- 
lations as a constraint among fields, we give a deriva- 
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tion of topological field theories which describe strongly 
interacting topological phases. Finally, we present our 



The pure gauge condition can be imposed by inserting 
the delta functional 



conclusions in Sec. VI 



II. FUNCTIONAL BOSONIZATION 

A. Functional bosonization 

We begin with a summary of functional bosonization. 
Our starting point is the fermionic partition function in 
D = d + \ space-time dimensions, 



n n 



(2.5) 



Z[A° 



V[^p,^]exp{iKF[^,^,A'^^]), (2.1) 



where Kp['i/},ip, A°^] is the fermionic action describing 
the (topological) insulator in question, and is an ex- 
ternal U(l) gauge field which couples minimally to the 
fermion field -i/jjip as a source (background). The corre- 
lation functions of the electrical current operator can be 
obtained from the generating functional as 



1 S 1 



in the path integral, where e^i' - a^ jg the Lcvi-Civita 
symbol in D space-time dimensions, and Il^t-'A^^ runs 
over D{D — l)/2 independent directions. For example, 
when D ~ 3 and D ~ 4, the delta functional is given 

by n. n'=i e'^-^'SifM), and 0. U'.T ^"""'"SifxAa]), 
respectively. The delta functional can be exponentiated 
by introducing an auxiliary rank (13 — 2) tensor field 
as 

Z[A''''] = J V[a,b]Z[A'''' + a] 

X exp ( - d^x-fo^,...e'^''-"^/„0[a]). (2.6) 

Using the invariance of the integration measure under a 
shift of the gauge fields, a -> a — A°^, we can write the 
partition function Z\A^^'^\ as 



I M-(xi) i M-(a;2) 



■lnZ[A<=^. (2.2) Z[A'''']= J V[a,b]Z[a] 



Functional bosonization makes use of the gauge invari- 
ance of the action and of the functional integral under a 
local U(l) gauge transformation: for a pure gauge field 
a, I, the partition function is invariant under the shift 

^ .4™ + a, 



X exp 



(2.7) 



Z[A'''' + a] = ZlA'' 
Thus, one can represent ZlA^^] as 



V[a]^^,,Z[A-^ - 



(2.3) 



(2.4) 



up to some normalization, where J 2? [a] pure is the path 
integral over the gauge field with the condition that 
it is a pure gauge. Thus, in terms of the field strength 
/pi/ [a] = d^a^ — d^a^, the allowed gauge field configura- 
tions are required to obey /^^ [a] = for all possible pairs 
of indices [fi < v). 

In the case of a system of fcrmions on an open man- 
ifold with fixed boundary conditions, such as a disk in 
two dimensions or a ball in three dimensions, this pro- 
cedure does not change the partition function. This is 
true independently of the dynamics of the fcrmions (pro- 
vided it is gauge-invariant). This is also true for a closed 
topologically trivial manifold (such as a sphere). On the 
other hand, in the case of fermions on a closed topo- 
logically non-trivial manifold, such as a torus in any di- 
mension, this procedure is equivalent to a averaging the 
partition function over the large gauge transformations 
on the torus. Thus, this is equivalent to averaging the 
partition function over "twisted" boundary conditions. 



From the partition function (2.7), the current correlation 
functions can be computed as the correlation functions 
of the tensor field h^y..., 

= (e'^i''i^--9,,6Ai...(a;i)e^^''=^--9,,6A,...(x2)---), 

(2.8) 

suggesting the correspondence (bosonization rule) 

f{x) ^ e^-^P-dMp-{x). (2.9) 

where the {D — 2)-form tensor field h^^... was introduced 
as a Lagrange multiplier field. This system has in fact 
a (local) gauge symmetry as its partition function is in- 
variant under 



(2.10) 



where the symbol {• • • } fully antisymmetrizes the indices. 
Observe also that in Eq. (^.7j ) the field strength for b^^... 
which is the {D — l)-form field h defined by h db (in 
the differential form notation), does not appear. 



A consequence of the bosonization rule of Eg. (2.9) and 
of the form of the partition function, Eq. (2^), is that 
local magnetic fluxes (or, in general dimension, magnetic 
flux "tubes" ) couple to the tensor field b. In addition, in 
a system with an energy gap the worldlines of fermionic 
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excitations (i.e. the quasiparticle currents) are minimally 
coupled to the gauge field a^. 

To complete the bosonization mapping, we need to 
evaluate the fermionic path integral Z[a\. Our discussion 
so far is quite general and applicable to gapped/gapless, 
interacting/non-interacting systems in any dimension. 
However, for a theory of massless fermions (free or in- 
teracting) the fermionic path integral Z[a\ can only be 
evaluated exactly in Z? = 1 + 1 dimensions. In general 
space-time dimension D the fermion path integral Z[a\ is 
a non-local (but gauge invariant) functional of the gauge 
field. On the other hand, if the band gap (mass) is fi- 
nite, an approximate form of the fermionic path integral 
Z[a\ can be obtained using the inverse mass expansion. 
Below, we will discuss first the case of non-interacting 
(topological) insulators and we will later extend these 
results to the interacting cases. The path integral cal- 
culation is most conveniently described by choosing the 
Dirac representative of topological insulators as a mi- 
croscopic model. The resulting bosonized action is then 
a sum of local gauge-invariant operators. In particular 
the parts of it which are topological, do not depend on 
microscopic details. These topological terms are always 
marginal and dominate the low-energy physics for d < 3, 
they also only depend on universal properties of the sys- 
tem. In particular they contain the information on the 
topological invariants of the microscopic model. 

In the next subsection we discuss some technical and 
subtle points about functional bosonization that might 
be skipped during a first reading. After this discussion 
we move on to discuss what this approach predicts in 
various dimensions D using several examples. 



B. Generalities 

The functional bosonization technique is rooted in 
the well-knpW|ft, fermi-bose equivalence in D = 1-1-1 
dimensions. Ejca (For reviews on bosonization see Refs. 
|32| , and Q). It is known that the density-density com- 
mutator of a Z? = 1 -|- 1-dimensional fermion is anomalous 
in the sense that it develops the so-called Schwinger term, 
while one would naively expect (for a system with a rel- 
ativistic spectrum) that the charge and current density 
operators at different positions (momenta) commute with 
each other at equal times. This is due to the underlying 
chiral anomaly. This anomalous commutator allows one 
to represent the fermion density operator in terms of a 
boson field. Indeed, if we represent the (normal-ordered) 
operators for the fermion charge density jo a-nd current 
density ji by a the two- vector = (j0iji)5 then the 
anomalous commutator [jo(a;), ji (y)] = {i/T:)dxS{x — y) 
implies that the current can be represented in terms of a 
scalar (Bose) field 0, i.e. = {l/y/n)€fj,i,d''(j), which is 
consistent with the requirement of local conservation of 
the charge current, = 0. In the Abelian case, the 

resulting effective field theory for the bose field <f> turns 
out to be free (up to a finite renormalization of the Lut- 



tinger parameter and of the speed of excitations due to 
interaction effects). In other words, the dynamics of the 
fermionic system (interacting or not) is fully represented 
in terms of the dynamics of its conserved currents. In this 
sense, in one space dimension, bosonization yields an ex- 
act hydrodynamic representation of the system. In fact, 
the full excitation spectra of both systems are identical 
and so are their partition functions. This exact fermi- 
bose equivalence is a consequence of the kinematic re- 
strictions of one-dimensional motion for systems with a 
relativistic (linear) dispersion. Using these identities one 
can establish a one-to-one mapping between seemingly 
different systems such as the D = 2-dimensional massive 
Thirringiptciodel and the D = 2-dimensional sine-Gordon 
modcl.EsH Such a correspondence also extends to Dj-^ 2- 
dimensional theories with non- Abelian symmetries.c3 

In 1-1-1 space-time dimensions these bosonization oper- 
ator identities can also be formulated in terms of the func- 
tional integral language. In this (functional) bosonization 
approach the fermion determinant of a Dirac operator 
coupled to gauge fields is computed by means of a local 
chiral transformation. At the classical level {i.e. in the 
Lagrangian) in 1+1 dimensions it is possible to (formally) 
cancel the coupling of the fermi field to a gauge field 
Afj^(x) by means of a local chiral transformation, which 
acts of the Dirac fermion as ipix) — >■ exp[i(l){x)j5]tp{x), 
where (f>{x) is an arbitrary smooth function of space-time 
coordinates x and 75 is a suitable Pauli matrix. In a 
gauge-invariant theory the chiral anomaly implies that 
the integration measure of the fermion path integral is 
not invariant under the chiral transformation. As a result 
there is a non-tpyial Jacobian associated with the chiral 
transformation.L^I This Jacobian is computable in terms 
of the fermion determinant which leads to the bosonized 
form of the effective action theory which is a local func- 
tion of the bose field (pix). This approach has been 
applied with great success to-aeveral 1 -f- 1-dimensional 
theories of massless fermions .cil In massive theories, the 
fermi-bose mapping relates the two descriptions as iden- 



3ation theory in 
3 It also yields an 



titles valid order by order in a p 
powers of the fermion mass term.L^ 
alternative derivation of non-AbelianJatOsonization yield- 
ing the Wess-Zumino-Witten modeled by meaes of an 
exact computation of the fermion determinant .e3 

For d > 1 it is not possible, in general, to find an ex- 
act mapping in the form of operator identities, between 
a local theory of fermions and a local theory of bosons. 
In practice, for d > 1 bosonization reduces to finding 
an effective field theory in terms of bose fields. Physi- 
cally this is largely due to the fact that the kinematic 
restrictions of one dimension do not exist for d > 1 and, 
as a result, the spectrum of a local fermionic theory, in- 
teracting or not, is iiot equivalent to the spectrum of a 
local bosonic theory.ES The resulting effective field the- 
ory is local only in the case of theories with a finite gap 
in the fermionic spectrum and for energies much smaller 
than the gap. The prime example of such a theory is the 
hydrodynamic field theory of the FQHE. 
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For the case of functional bosonization, the problem 
is that the resulting fermion determinant for massless 
fermions coupled to gauge fields is a gauge- invariant but 
non-local functional of the gauge field and, contrary to 
the case of D = 2, cannot be computed in closed form. 
However in the case of massive fermions the contribu- 
tion of the fermion determinant to the effective action 
of the gauge fields can be expressed in terms of local 
gauge-invariant terms (for manifolds without boundary) 
with coefficients that involve powers of the inverse of the 
mass, i.e. a gradient expansion. This approach natu- 
rally is only meaningful for massive theories. In the case 
of massive Dirac, fermions in D = 2 -I- 1 dimensions it 
has been shownEltj that the correlation functions of the 
conserved fermionic currents is the same as the correla- 
tion functions of a dual topological Chcrn-Simons gauge 
theory in the low energy regime (only!). A mapping of 
the correlators of other fermion bilinears can also be de- 
termined but, again, only in the low-energy and long dis- 
tance regime. 

Below, we will show that functional bosonization can 
be used to derive hydrodynamic theories of topological 
insulators in d = 1, 2, 3, 4. Specifically, we give the func- 
tional bosonization results for non-interacting topological 
insulators in D = d + 1 with d = 1 , . . . , 4 which crucially 
conserve electromagnetic U(l) charge. These topological 
phases belongs to "the complex symmetry classes" and 
to "the primary scries" of the eight real symmetry classes 
in the periodic table; By "complex symmetry classes" we 
refer to symmetry classes A and AIII, while the "pri- 
mary series" of topological insulators (superconductors) 
is located on the diagonal of the periodic table, and are 
characterized by an integer Z topological invariant. 

After this we will discuss field theory descriptions of 
non-interacting topological phases characterized by Z2 
invariants. These theories can be obtained from the 
"the primary series" by dimensional reduction and can 
be characterized in two different classes: first and sec- 
ond descendants. These phases include the time-reversal 
symmetric topological insulators in £) = 3 -I- 1 (a first 
descendant) and D = 2 -|- 1 (the quantum spin Hall ef- 
fect, a second descendant) dimensions. It is also possible 
to extend our discussion to topological phases (topolog- 
ical superconductors) which cousckm; non-Abelian cur- 
rents such as a spin SU(2) current.E3 Such non-Abelian 
functional bosonization is discussed in Appendix ^ 



Altland-Zirnbauer classification,!^ they belong to sym- 
metry class AIII. In D = 1-1-1, topological insulators 
in symmetry class AIII are characterized by an integer- 
valued topological invariant, the winding number J^i— lA 
canonical example for such systems is polyacetylene.E3 

In addition to the sublattice symmetry, we can further 
impose time-reversal symmetry which squares to be ei- 
ther +1 or —1 (corresponding to symmetry class BDI or 
CII). In cither case, band insulators in _D = 1 + 1 dimen- 
sions in these symmetry classes are characterized by an 
integral topological invariant v, in much the same way 
as symmetry class AIII in = 1 + 1. Below, we will 
consider symmetry class BDI in Z? = 1 + 1, which is in 
the primary series in the periodic table. 

The symmetry class BDI can be realized eithep|i|t-tpS|ms 
of complex (Dirac) or real (Majorana) fermions. OEEl In 
the latter case, instead of imposing a sublattice symme- 
try, we would impose a reality condition (originating from 
the fact that we are dealing with Majorana fermions) 
combined with time-reversal symmetry. The reality con- 
dition is equivalent to a charge-conjugation symmetry 
and when combined with time-reversal plays a role sim- 
ilar to sublattice symmetry. Our discussion below, how- 
ever, will focus on the case of complex fermions, as func- 
tional bosonization takes advantage of the presence of a 
U(l) global symmetry. 

To summarize, wc consider (topological) band insula- 
tors with U(l) global symmetry (particle number con- 
servation) in symmetry class AIII or BDI i n D = 1 + 1 
dimensions. For these cases, Z[a] in Eq. (2^) can be 
evaluated as (in the Euclidean signature) 



lnZ[a] = 



i0 
2^ 



d^xe^'^dna 



(2.11) 



Here, the angle 9 is given in terms of the bulk topological 
invariant, the winding number v, as 9 = vk mod 27r. In 
Eq. ( 2.11 ), terms with more fields and derivatives such as 
Maxwell terms are suppr essed as they are more irrelevant 
[indicated by • • • in Eq. (2.11)]. 

Thus, the bosonized partition function is given by 

= j 2?[a,5]expi j d^xC (2.12) 
with the effective low-energy Lagrangian density 



C. Functional Bosonization Examples 

1. D = 1 + 1 (AIII or BDI) 

There are no non-trivial topological insulator states 
in = 1 + 1 dimensions if we do not require a pro- 
tected symmetry. That is, topological insulators (in fact, 
topological phases in general) in _D = 1 + 1 must be 
symmetry-protected ones. In the following, we consider 
band insulators with sublattice (chiral) symmetry. In the 



C 



(2.13) 



In summary, the effective low-energy theory of massive 
fermions in 1 + 1 dimensions is a BE theory (the first 
term) and a topological term whose coupling constant is 
the topological invariant 9 = ■nv where u is the winding 
number. Notice that \n D = 1 + 1 dimensions the BE 
theory involves a vector (gauge) field and a scalar 
field b. 
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2. D ^2+1 (A or D) 

In Z) = 2 + 1, we discuss the bosonization of topological 
insulators in symmetry classes A and D. They are topo- 
logical insulators belonging to the primary series in the 
terminology introduced above. These topological insula- 
tors are characterized by an integer topological invariant, 
Ch, the first Chern number, which is nothing but the Hall 
conductivity a^y (in units of e^/h). The symmetry class 
A is defined as a set of fermionic quadratic Hamiltonians 
which possess no symmetry, and the canonical example 
of the topological insulator in this class is the integer 
QHE. The symmetry class D is obtained, from symme- 
try class A, by imposing additional particle-hole symrr^ 
try while keeping the electromagnetic U(l) symmetry.Eil 
This should be distinguished from the superconducting 
realization of symmetry class D, where there is no con- 
served U(l) charge whereas the particle-hole "symmetry" 
simply reflects the fact that fermionic Bogoliubov quasi- 
particles are a real fermion. 

For these cases, Z[a\ can be computed in the low- 
energy limit (again in the Euclidean signature) as 

In Z[a] = ^ y d^x e^^'^a^a^flA + • • • , (2.14) 

where Ch is the integer- valued bulk topological invariant, 
the Chcrn-numbcr, of the topological insulator. As in the 
case oi D = 2, terms which are less relevant compared 
the Chern-Simons term, including the Maxwell term, are 
suppressed by powers of the inverse of the mass (the en- 
ergy gap). Thus, 

= J V[a,b]cxpi J d^xC, (2.15) 

with the effective low-energy Lagrangian density 

Ch 



C 



AT) + —e>'''\,,d,ax. (2.16) 
47r 



Thus, here too we obtain a BF theory (this time in 2 + 1 
dimensions) and a topological invariant term, the Chern- 
Simons action, whose coupling constant (or "level") is 
the (topological invariant) Chern number of the micro- 
scopic band. In D = 2-1-1 the BF term relates two 
gauge fields, &^ and a^. This effective action has a formal 
similarity with the hydrodynamic theory of the FQHE. 
However in the hydrodynamic theory of the FQHE the 
two gauge fields are the hydrodynamic field 6^, and the 
statistical gauge field a^, with the important difference 
that Chern-Simons term affects the hydrodynamic 
ficld.r'Trn In the present case the bosonized action, Eq. 



( 2.16 ), represents an integer QHE with a quantized Hall 



conductance cixy = {e^/h)Ch. Thus, this is a theory of 
the quantized anomalous Hall effect. 

3. = 3-1-1 (AIII or Dili classes) 

In D = 3 + 1, topological insulators in symmetry 
class AIII are characterized by an integral topological 



invariant, the three-dimensional integer-valued winding 
number and protected by a chiral symmetry. A lat- 
tice tight-binding example of topological insulators in 
D = 3 + 1 in symmetry class AIII can be found in Ref. 
|6^ . This topological insulator is somewhat analogous to 
the time-reversal symmetric topological insulator in sym- 
metry class AH, in that it supports a Dirac fermion sur- 
face state, and has non-trivial axion response to the ex- 
ternal electromagnetic field. The difference is, however, 
that the latter is characterized by a Z2 topological in- 
variant, rather than an integer topological invariant. We 
will discuss the bosonization of the time-reversal sym- 
metric topological insulator in symmetry class AH in a 
later section. 

The topological insulator in symmetry class Dili is 
similar to the topological insulator in symmetry class 
AIII above but with the requirement of an additional 
particle-hole or time-reversal symmetry (which lead to 
equivalent results since these symmetries, when combined 
with the sub lattice symmetry give the third). As in the 
cases of D = 2 and D = 3 above, this realization of sym- 
metry class Dili in terms of complex fermions {i.e., with 
conserved U(l) charge) should be distinguished from the 
superconducting realization of symmetry class DHL 

For these cases, Z[a] can be computed, again in the 
low-energy limit, as (in the Euclidean signature) (see, for 
example, Ref. ^ for calculations in terms of the Dirac 
representative) , 

\nZ[a]^^J d^x t^^^Pd^^a^dxap + • • • . (2.17) 

Here the angle ^ is related to the winding number as 
d = v'K mod ItxVA Thus, 



Z[yl'= 



X>[a, 6]expi / d^ xH, 



(2.18) 



with the effective Lagrangian 



9^ 



jd^a.d^a" + ■■■ 



(2.19) 



In the last line, we have written down the Maxwell term 
explicitly since it is also marginal in I? = 4. 

4. D = A + l (A or All classes) 

Finally, we discuss the bosonization of the topologi- 
cal insulators in D = 4 + 1. In fact, the bosonization 
rule applies to any dimension. One of our motivations 
to study the case of D = 4 + 1 is that the lower di- 
mensional topological insulators, in particular, the time- 
reversal symmetric topological insulators in Z? = 3 + 1 
and D = 2 + 1, are the first and second desccndaiiis of 
the _D = 4 + 1 topological insulator in class AILO-E-a In 
D — 4 + 1, topological insulators in symmetry class A 
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and All are characterized by an integer-valued topologi- 
cal invariant, the second Chern number Ch2. 

For these cases, Z[a] can once again be computed as 
(in the Euclidean signature) 



\nZ[a] ='^J d^xe^^-'^P-a^d.a^dpa, 



1 



du,aud^a!' + • • • 



(2.20) 



where the integer Ch2 is the second Chern number of the 
topological insulator. Thus, 

= j 2?[a,6]expi j d^xL (2.21) 

where the effective Lagrangian is now given by 



(2.22) 



where we have dropped the Maxwell term in the last 
equation. 



5. Interactions 

We have not included interactions so far, but pertur- 
bative effects of weak interactions can easily be taken 
into account. As an example, consider the local current- 
current interaction (which is the analog of the Luttingcr- 
Thirring interaction in _D = 1 + 1 dimensions) 



(2.23) 



where u is the coupling constant. Because of the 
bosonization rule, this leads to a term of the form 

-u j (i^x(e^'^^^ -9.6Ap...)(eMo/3-5"&^^ -)- (2-24) 

That is, the interactions are represented by a Maxwell- 
like-term for the h field. 



III. PHYSICAL PICTURES FROM 
FUNCTIONAL BOSONIZATION 

In this section, we discuss the physical picture that 
functional bosonization provides in three different situa- 
tions with increasing complexity: an example of a topo- 
logically trivial insulator, a time-reversal breaking topo- 
logical insulator (the QHE or quantum anomalous Hall 
effect) in D = 2 -|- 1 dimensions, and a time-reversal in- 
variant topological insulator in Z? = 3 + 1 dimensions. 



A. Topologically trivial insulators 

We start our discussion from topologically trivial in- 
sulators and we will make a few comments. We will dis- 
cuss: (i) the normalization of the gauge fields when the 
BF theory description is considered on a compact spatial 
manifold (ii) functional bosonization and the dual picture 
of insulators, and (iii) the presence/absence of fermionic 
excitations in the low-energy spectrum. In fact, these 
comments are not limited to topologically trivial insu- 
lators, and we will continue our discussion of them for 
topological insulators in the later sections. 

(i) Functional bosonization gives us a bosonized de- 
scription of non-topological band insulators in terms of 
the BF theory without a CS or 9 term, i.e., the pure 
BF theory. In = 2 -|- 1 dimensions, the bosonized La- 
grangian of a non-topological band insulator is given by 
Eq. (2.16) with vanishing Chern number, Ch = 0. On 
a compact spatial manifold, such as a torus, we should 
be careful of the normalization of the gauge fields which 
must be done relative to their compactification condition. 
On the torus, where the spatial coordinates {xi^X2) are 
identified periodically + L.i, the gauge fields a) 



("i = bf, and al 



i), in the ctg = gauge, are normal- 



ized such that they are global gauge-transform equivalent 
to + 2Tmj / Lj where n-j=i,2 are integers. When prop- 
erly normalized, the BF Lagrangian is given by 



2k 



An 



(3.1) 



with k — 1. With this normalization, the ground state 
degeneracy is one {i.e. no ground state degeneracy). 

(ii) The BF theory description of (non-topological) 
band insulators psjembles the dual description of BCS 
superconductors. The bosonized theory of the 
trivial-band insulator from functional bosonization is 
naturally pelatcd to the dual theory of (band/Mott) 
insulators. EllGa Below, let us take a closer look at this 
by focusing on _D = 2 + 1, although a similar discussion 
applies for _D 7^ 2 + 1. 

In (2-|-l)d superconductors, a defect (vortex) of the 
condensate A(r) traps a magnetic field V x Aem(^); 
the phase 4>{f) of the superconducting order parameter 
A(r) = |A(r^| expi(/)(r) is related to the magnetic field as 
■^<^o §gg\J4'-dl = /g(VxAEM)-c'S', where <I>o is the flux 
quantum, and if we reinstate ft., c, and e, <f>o = ^c/(2e); 
the magnetic flux is localized in the region S. 

In the dual description of an insulator, we view the 
charge density modulation 5p[r) caused by external dop- 
ing as a point-like "defect". As in the (2-|-l)d supercon- 
ductor, we postulate that such a defect traps a flctitious 
magnetic field V x a 

( SpdS I {\7 xa)-dS, (3.2) 
Js Js 

where is a fictitious gauge field. As it will become 
clear, this gauge field is the same gauge field appearing 
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in Eq. (3.1). As the local charge density modulation is 
viewed as a vortex, we can postulate the existence of an 
order parameter $. The phase of <i>, $ = |$|cxpi8, is 
related to the charge density as 



ve • dl 



dS 



SpdS. 



(3.3) 



As in the dual theory of superconductors, the above 
description in terms and $ can be dualized (in the 
condensed phase) and can instead be rewritten in terms 
of two gauge fields and 6^ , which are the gauge fields 



appearing in Eq. {'3A). Upon dualization, 6^ couples 



minimally to the vortex current of $, jy = e^'^^di.dxQ. 
This is in fact consistent with the bosonization rule 

(iii) The bosonic descriptions (in terms of and $, or 
in terms of and 6^) do not contain any fcrmionic exci- 
tations even though we started from a gapped fermionic 
theory (a band insulator). Drawing again an analogy to 
BCS superconductors, this resembles to the fact that, un- 
der duality, a superconductor in a uniform field is equiv- 
alent to the insulating phase pf JJae same bosonic the- 
ory at uniform charge density.EjO In the latter theory, 
the charges are gapped Cooper pairs, not fermions. For 
this reason the bosonic theory only describes arrays of 
Josephson junctions in which the fermionic excitations 
do no exist (they are bound on each superconducting 
grain). The connection between the gapped fermionic 
theory and the bosonized descriptions is that they share 
the same topological limit in the form of a BF theory. As 
we will see in the next section, the bosonized description 
of a D = 2 -|- 1 dimensional topological insulator is dif- 
ferent and supports fermionic excitations (electrons) due 
to the presence of the Chern-Simons term. 



a Chern insulator. In fact, this effective field theory de- 
scribes a topological fluid. This can be checked by noting 
that it describes a system with a non-trivial ground state 
degeneracy k^, where g is the number of handles of the 
surface. In other terms, a theory with k 7^ 1 describes 
a fractionalized topological insulator, e.g. a fractional 
quantum Hall fluid, as w e will see in Section ^ 

The Lagrangian (3.4) is the hydrodynamic effective 
field theory for _D 2 -I- 1 dimensional topological insula- 
tors. Some examples include the QHE realized in a two- 
dimensional electron gas with Landau levels produced by 
a uniform magnetic field, as well as the Chern-insulator 
i.e., lattice fermion systems with non-zero Chern number 
but without uniform magnetic field, such as t he H aldane 
modcO). Note that we have arrived at Eq. ( 2.16 ) with- 
out using composite particle theories (see below for more 
comments); and that for the Chern- insulator, the flux 
attachment transformation is highly non-trivial but not 
impossible. 

It is readily seen that the bosonized theory repro- 
duces the QHE: From the equations of motion, 5S/5a^ — 
5S/5h^ = 0, 



— e*" dt,bx + —e^ a^ax 
in zn 

(ax - A^^) = 0. 



0, 



From the dcflnition of the electrical current 
SS 2k 



in 



we then conclude the QHE as 



27r 



(3.6) 



(3.7) 



(3.8) 



B. Topological insulator in D = 2 + 1 

Let us now move on to a more non-trivial discussion of 
a topological insulator. The bosonized Lagrangian of the 
D = 2+1 dimensional topologi cal in sulator (in symmetry 
class A or D) is given by Eq. ( 2.16 ): 



2k 



Ch 



C = -—b^e''''^d,{ax - AT) + —e^-^a^,d,ax 



in 



in 



= + ^h^^^'^^d^AT, (3.4) 

where we integrated by parts once, and introduced the 
following notation: (a^,Q!^) — {bf^^a^^) and 



K = 



-k 
-k Ch 



(3.5) 



We have normalized 6^ as before (assuming we are on a 
compact manifold such as a torus) and k = 1 . 

This effective theory also makes sense for k 7^ 1. How- 
ever in this case this effective hydrodynamic theory can- 
not be derived from a free-fermion system, even if it is 



While functional bosonization delivers the two- 
component CS theory (the BF-CS theory), it is possi- 
ble to derive the familiar single-component CS theory by 
integrating over a^: If we eliminate by using the equa- 
tion of motion, for the simplest case of k = Ch = 1, we 
obtain the effective action 



in Zn 



(3.9) 



This is the familiar single-component CS theory (here for 
the case of the integer QHE). Observe that, as it should, 
the CS coefficient in front of e^^^b^j^d^bx has the oppo- 
site sign as compared to the CS theory of the response, 
(Ch/47r)e^''^A"^a^^^''. 



The transition from the BF-Chern Sim ons-theory (3.4) 
to the single-component CS theory (3.9) is quite analo- 
gous to the derivation of the CS description for the (in- 
teger) QHE in terms of the composite boson theory. In 
the composite boson theory, one introduces a CS the- 
ory which attaches a flux (fluxes) to electrons, and the 
bosonic composite particle (electron -I- flux), in terms 
of a dualized language, is expressed by a U(l) gauge 
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field; These two U(l) gauge fields correspond to and 
&^ in functional bosonization. In fact, if we "dualize 



back" the BF-Chern Simons-theory (3.4), we obtain the 
Chern-Siomns-Landau-Ginzburg action, in terms of the 
CS gauge field and and some bosonic field $. Com- 
paring with the composite boson theory, the bosonic $ 
could be interpreted as a composite boson field. Thus, in 
functional bosonization, while we have not used the fiux 
attachment, the resulting effective field theory is quite 
similar to the composite boson theory of the IQHE. 

As clear from the comparison to the composite boson 
theory, the hydrodynamic theory ( ^.9[ ) [and hence ( ^.4[ )] 
includes a fermionic excitation (electron) due to the pres- 
ence of the Chern-Simons term. The world-lines of point 
defects (point sources for the vortex current that couples 
to 6^), when linked, pick up a tt phase (i.e., a fermionic 
sign) due to the the Chern-Simons term. This should 
be contrasted with the hy drod ynamic theory for topo- 
logically trivial insulators ( ^.l[ ) which does not have the 
Chern-Simons term, and thus no fermion statistics. 



C. Topological insulator in D = 3 -I- 1 

The bosonized Lagrangian of the D = 3 + 1 dimen- 
sional (topological) band insulator in symmetry class 
AIII (Dili) is given by Eq. ( ^.191 ): 



C 



87r2 



(3.10) 



where we have dropped the Maxwell term for simplicity. 
The theta angle 9 is i^tt where i/ G Z is the topological in- 
variant (winding number). In the following, we will con- 
sider a background 0-field which varies in space (and time 
when necessary), 9 = 9{x). This allows, for example, 
an interface between band insulators with different topo- 
logical charges ly, or a. topological band insulator where 
symmetry conditions defining a given Altland-Zirnbauer 
symmetry class (such as sublattice, particle-hole, time- 
reversal symmetries) are (weakly) broken in the bulk. 

The effective action ( 3.10| ) obtained from functional 
bosonization was discussed previously in Refs. and 
|7^ . The effective action ( 3.10|) should be compared with 
the one proposed in Ref. 24|: 

Zn Zn 



(3.11) 



where C ~ ±l/(87r). One can check, upon integration 
over a and b, we reproduce the axion term 



(3.12) 



As compared to Eq. (3.10) the effective action (3.11) does 
not have the axion term for field. However, assuming 



there is no monopole in fie ld co nfigurations, with the 
gauge transformation in Eq. ( ^.10 ) 



IBtt' 



(3.13) 



one can transform Eq. ( |3.10| ) into Eq. ( |3.1l| ). 

As in the case of D = 2 -I- 1, one can readily verify that 
the bosonized effective Lagrangian reproduces physics of 
the topological insulator m D = 3 + 1 dimensions: By 
eliminating b^i, and by making use of equations of 
motion, SS/Sa^ = SS/Sb^i, — 0, 

e^^^'dMp - -^^''"^'du {9dxap) = 0, 
e^''^^dx{a,-Al^)^0, (3.14) 
one obtains the axion term for the external gauge field 
1 



87r2 



(3.15) 



as expected for the response of the topological insulator 
in symmetry classes AIII and Dili in D = 4. Also, from 
the bosonization rule = e^^^^dybxp, the electrical cur- 
rent is given by 



r 



ID 



(3.16) 



Taking to be a step function which jumps from tt — > 
at the system's boundary, this equation tells us that 
surface Hall conductance is am = t- — ttt where we have 

■^y 47r 2h 

reinserted units. Alternatively, when 9 = t: = const., but 
has a monopole configuration, this equation tells us 
the monopole acquiijeSja. charge which is a manifestation 
of the Witten effect Jia'ta 

To discuss the effect of monopoles, it is convenient 
to introduce a monopole gaug e field f^y[A'^^] and /^[a] 
for A°^ and a^, respectively.EZlOlZj In the presence of a 
monopole in A'^^, the monopole gauge field .f^y[A°^] is a 
field which is related to the monopole current j^^'^ as 



2 " " IX J Ap L 

More specifically, it is given by 
1 



(3.17) 



CW^'Mx) = -e,,xp5^P{x- S) = ~8p,{x- S), (3.18) 
where S\p{x; S) is singular on the world-surface S, 



S^'^ix-^S) 



dcrdr 

X S^^^lx^ 



(9X^(cr, r) dX'^ia, r) 



da 
X{a,r)] 



dr 



(3.19) 



The surface S is the world-sheet of the Dirac string, and 
which is bounded by a world-line L of the monopole, 
dS = L. One verifies 

le^^''^f^Jxp{x;S)=SP{x■,L), (3.20) 
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where 



da ^^y^^h ^^-^ix-Xia)]. (3.21) 

OCT 



One then conchidcs /^[^™] is related to the monopole 
current as Eq. ( ^.17 ). 

To describe the presence of a monopole in terms of the 
monopole gauge field, we separate the physically observ- 
able field strength in the Lagrangian into the monopole 
free part and the monopole gauge field, ffj,^[a] — > /pi/[a] — 



327r2 



(/apH-/apN) 



(3.22) 



The Lagrangian is invariant under the ordinary U(l) 
gauge transformation. 



Jp ^ + 9pA. 



(3.23) 



In addition, it is also invariant under the monopole gauge 
transformation, 



C[a]^Ci<^]+d^V.~d.r,^, 



af_c ^ afj_ + rj^. 



(3.24) 



This second gauge invariance originates from the fact 
that for a given Dirac string, there are many surfaces 
satisfying dS = L. 

Let us now focus on the case where 9 is constant 
throughout the bulk. We can use the equations of motion 
SS/Sofj, = SS/Sbfj^i, = derived from the effective action 
( |3.22 ) to find that the electrical current is given in terms 
of the monopole current as 



(3.25) 



Given that in D = 2 + 1 the gauge field plays the 
role of attaching a flux (fluxes) to electrons, it is tempt- 
ing, and perhaps instructive, to view the gauge field 
flp in _D = 3 + 1 as a proper generalization of the CS 
"statistical gauge field." From the equation of motion 
e'"'^Pd^b),p = ~{0/4■K^)j^'^', the gauge field is attach- 
ing a monopole to the electron (and as in = 2 + 1, 
the temporal component ao enforces a constraint). If 
we further choose, in the presence of a vortex line, the 
worldsheet of the Dirac string to be identical to the vor- 
tex worldsheet, the boundary of the vortex worldsheet 
(the end of the vortex) is the worldline of a mo nopole 
{i.e., a dyon since 9 ^ 0). In this case, Eq. (B.2E) gives 



h^,^ = -(6l/87r2)/M [A<=''] up to a choice of gauge. 



IV. DIMENSIONAL REDUCTION 

We have derived the effective BF topological field theo- 
ries for the primary series of the topological insulators in 
the periodic table. In microscopic fermionic theories, TL^ 
topological insulators (such as the time-reversal symmet- 
ric Z2 topological insulator in Z? = 3 -I- 1 dimensions, and 
the quantum spin Hall effect in D = 2-1-1 dimensions) can 
be derived as a "descendant" of the pciHiary series by the 
Kaluza-Klein dimensional reduction. Ej We now discuss 
the Kaluza-Klein dimensional reduction of the BF type 
effective action we derived. This is expected to deliver an 
effective bosonic field theory description for descendant 
Z2 topological insulators. Similar dimensional reduction 
was discussed in Ref. |l^ for the topological response the- 
ories. 

Let us start from the D = 4 + 1 dimensional BF-CS 
theory (the "parent" theory): 



Ch, 



247r2 



d''xe^"'^P''af,d,axdpaa + 



(4.1) 



where pL,v,\, 
as 



0, 1, . . . , 4. We now expand the fields 



+00 



<^{xi,n^) (4.2) 



n^,— — 00 



where <f> represents a field, Ly^ is the circumference of 
direction, and Xi = {xo,xi,X2,x^) and x^ = w. By 
"shrinking" the direction by taking — > 0, the "par- 
ent" D = 4 + 1 dimensional theory is reduced to a descen- 
dant theory in Z? = 3+1 dimensions. In doing so, we view 
each Fourier mode ^{xi, n„) as a field in Z) = 3 + 1 di- 
mensions (Kaluza-Klein modes). The modes with ^ 
have a gap which grows as we take L.^ — >■ and hence 
at low energies only the Fourier modes with ~ are 
important. [This can be seen in the presence of proper 
kinetic terms (e.g., the Maxwell term) for and b^^x, 
which are not shown explicitly above] . We define 

a^,{xi,n^ = 0) =: (j){xi)/L^, 



AZ{x„n^,=Q) =:9-^{x,)/L^, 



(4.3) 



{i,j = 0,...,4). The resulting Z? = 3 + 1-dimensional 
Lagrangian inherits the field content of the parent theory: 
one form gauge fields ai and Af^, a three form gauge field 
6ijfc, where the space-time index in Z? = 4 dimensions 
i,j,k runs from to 4. In addition, as introduced in 
Eq. (4.3), it also contains two scalar fields, </> and 9'^^, 



and one 2-form gauge field Uij. The scalar field 9'^'^ is 
the background axion field, which satisfies 9'^^ = 7r(0) 
inside (outside) of the Z3 = 3 + 1 dimensional topological 
insulator. Putting together, the effective topological field 
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theory for the D = 3+1 dimensional topological insulator 
is given by the Lagrangian, 



2-K 



27r 



(4.4) 



D = 3 + 1 dimensional fields, as 



bzij{Xi,Tlz) ='■ '^gij{xi) / Lz, 

Uzi{xi,nz) =: ^Vi{xi)/Lz, 



(4.7) 



where the space-time index ^, v, A now runs from to 4. 
The last two terms resemble the effective field theory of 
the D = 3 + 1 dimensional topological insulator in the 
prim ary s cries, i.e.,, the BF theory with the axion term , 
Eq. ( ^.ISP , whereas the first term is absent in Eq. ( 2.19 ). 
The Z2 nature of the system lies in the restriction on 
On20°^ which is quantized by time-reversal symmetry to 
be fixed at Ch26'°'' = 2mi or (2n -I- l)7r where n is an 
integer. The Z2 nature is manifest in the fact that the 
bulk of the material only uniquely determines whether 
0^28'^^ is. an even or odd multiple of tt not what the value 
of n is.U If wc calculate the electrical current response 
we find 



Sb^Xp 2tt 



(4,5) 



which gives rise to 



Ch2 

47r2" 



-dxA'^;' 



(4.6) 



where we have ignored possible monopole contributions 
to the current. This response implies that if Ch29°^ = 
(2n + 1)tt inside the material and 2mT outside then there 
is a half-integer quantum Hall effect on the surface. 

We can continue the reduction down to £> = 2 + 1 to 
obtain a hydrodynamic effective field theory for the quan- 
tum spin Hall effect; we separate the D ~ 4-dimensional 
coordinates (ccq, xi, a;2, aJa) into the D ~ 3-dimensional 
ones (a;o,a;i,X2) and a X3 = z. The z direction is com- 
pactified on a circle with radius and then we take 
Lz — >■ 0. The fields can be Fourier decomposed in the 
^-direction, as in Eq. (O), and only the Fourier modes 
with Uz = (the momentum quantum number in the 
z-direction) are kept. This second step of dimensional 
reduction introduces, in addition to two new scalar fields 
ipjX"^^, a vector field Vi, and a second 2-form gauge field 
9ij = 0,1,2,3). They arc defined, from of the 



The resulting D = 2 + 1 dimensional Lagrangian is 

-(a^ - Al^)d,vx + ^e^'^VS^V'S.aA, 



2TT ' " > 



(4.8) 



where the space-time index /i, i/, A now runs from to 3. 
We recognize the third term £^'^'''(0^1 — A™)9,yUA as the 
BF coupling in _D = 2 + 1 dimensions with 6^ replaced by 
Vfj_ . The presence of such term is largely expected based 
upon the phenomenology of the non-chiral (helical) edge 
modes of the quantum spin Hall effect. However, there 
are some additional terms which have not been previously 
discussed in dynamical gauge theories: the coupling to 
the real scalar fields and ip which could be combined 
into a single complex scalar field. We note that this the- 
ory, unlike othot byjdrodynamic theories of the quantum 
spin Hall effectErEJ only has a single U(l) gauge invari- 
ance (unlike the more conventional U(l) x U(l) for charge 
and spin). This is a natural result since in real materi- 
als spin is not conserved generically and only the local 
charge U(l) invariance is preserved. 

It is a simple exercise to read off, from the above effec- 
tive field theories, a response of the system to the external 
field. For example, for D = 2 -I- 1 dimensions, the def- 
inition of the electrical current j'^ := SS/SA^^ together 
with th e eq uation of motions derived from the effective 
action (^^sj) , 



5S 

5S 
Sa,, 



gives rise to 



p.uX 

= .f^—d.vx, 



2tt 



~2^ 



Ch2^ 



^--^a. (.^^A^) = 0, (4.9) 



(4.10) 



If we assume that e^^dp,dv'4' = then we can write this 
in terms of the external scalar fields as 



f = ^e^^^^d^e'^^dxx' 

477'' 



(4.11) 



A similar expression for the electrical current was derived 
in Rcf. ^ in terms of the effective topological response 
theory. If we think of the field 4>dxi} as a velocity field 
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V\ then this response imphes that there is charge bound 
to the vorticity of the velocity field j° = f^V x V, sim- 
ilar to the quantum Hall effect where charge is bound to 
the magnetic flux. A simple example is the case when 
6»<==^(x,y) = ne{y) and X {x,y) = 2n{e{x) - 1/2) where 
9'^^ represents a jump in the the axion angle at an edge 
andx'^'^ represents a magnetic domain wall on the edge. 
The charge confined to the edge magnetic domain wall 
is Q ^ / dxdy27r^S{x)Siy) = The distinction 

between integer and half-integer charge, i.e. Ch2 even or 
odd, gives an electromagnetic characteristic to determine 
the Z2 topological naturcj-of the 2d time-reversal invari- 
ant topological insulator.M Another consequence of the 
action is the fermionic mutual statistics between quasi- 
particles that couple as j'^a^ and i^^w^. When Ch2 van- 
ishes the currents j^, have mutual fermionic statistics 
but when Ch2 7^ the current is shifted such that 
and K^^ carry mutual fermionic statistics where 

~f^f-^e^^''''^,{<i,^^i,). (4.12) 



V. FRACTIONAL STATES 

Our discussions so far concern a bosonized description 
of non-interacting topological band insulators, or weakly 
interacting topological insulators that are adiabatically 
connected to a topological band insulator. In these sys- 
tems, functional bosonization gives rise to the BF topo- 
logical field theory with the unit level k = 1. In this sec- 
tion, we will discuss the possibility of topological states 
that can be induced by strong interactions and charac- 
terized by, e.g., non-zero ground state degeneracy. While 
the existence of strongly interacting topological insula- 
tors with topological order is not well established mi- 
croscopically, we can use the hydrodynamic bosonization 
formalism of the preceding sections to explore possible 
mechanisms that lead to a topological order. One such 
route is the fractionalization of electrons due to strong 
correlations. We will implement tbisj-through the par- 
ton construction a la Blok and Wen,H&3 to find effective 
topological field theories. 



Parton Construction for the Chern Insulator in 

D = 2 + l 



statistics. Since the electrons do not split in reality, we 
impose an constraint 



]^=3\i\ *,j = i,...,p+i. 



(5.1) 



on the parton densities jfi . The i-th parton carries elec- 
tric charge e^, and e = = 1- We assume the all 
partons are in their integer quantum Hall states, but we 
treat the i = 1, . . . , p-th flavors and the i ~ p + 1-st flavor 
differently. For i = 1, . . . , p, we assume the flUing faction 
is v'^'^^ = 1 whereas for i = p + 1, z^*^'-* = m. The total 
filling fraction is 



1 



m 



p + 1/m mp + 1 



(5.2) 



Each parton can be bosonized by functional bosoniza- 
tion. For i — 1, . . . , p-th partons, they can be described 
by the following D = 2 -|- 1 BF theories with the CS term, 

= ^e6(*)0a(*) + —ea'^'^da'^'^ - ed^'^ ■ A'''', (5.3) 
2tt 47r tj , y J 

where the repeated indices i are not summed, and we 
have introduced a short hand notations ■ j = 
and eada = e^^^^a^d^a\, etc. On the other hand, for 
the i = p + 1-st parton, since we have m filled Landau 
levels, we introduce m separate gauge fields, and 



a 



1=1,. ..,m 



each representing the condensate in the /-th 
Landau level. The Lagrangian for the i = p -|- 1-st parton 
is 



£(P+1) ^ ^ 



1=1 



-1 

2^' 



1 

47r' 



(5.4) 



The total Lagrangian is given by £ = Yl\=i ■ 
With the constraints imposed on the parton densities, 

e^'^'^d^bf = e'^^^9^6^P+^' = ^ t^^^d^bi, (5.5) 



{i = l,...,p), the total Lagrangian (after solving the 
constraints) is given by 



0) 



27r ^ 27r ^ 

47r ^ 47r ^ 



(5.6) 



We start with a construction of time-reversal break- 
ing fractional states in Z? = 2-1-1 dimensions by 
combining functional bosonization and the parton con- 
struction. This approach is equivalent to the compos- 
ite parti cki / pomi ysitc boson and composite fcrmion) 
theories. E§E§LlrE3L3 We focus on the effective field the- 
ory for the hierarchy states at v = m/(mp-|- 1) in terms 
of the parton construction. Following Refs. ^ and ^ 
we first split the electron into p -I- 1 partons. Here, p is 
an even integer, and we require partons to obey Fermi 



(We have switched off A™ for simplicity.) This can be 
written more compactly as 



47r 



(5.7) 



where a = {b^ , a^*)), and the A'-matrix is 


















(5 


~ p,m 





t 
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where is an to x to identity matrix and Jm.p is the 
m X p matrix with all matrix elements = 1. 

The topological order encoded in the JiT-matrix CS 
theory (5.7) can equivalently be described by a multi- 
component CS theory with fewer components. For ex- 
ample, when m = 1 and p = 2, the filling fraction is 
1/ = 1/3 and 



K = 




(5.9) 



The ground state degeneracy can be read off from 
|Det.ftr| = 3, as expected for the Laughhn state at 
v = 1/3. We can integrate a*^*^ one by one to ar- 
rive at the more familiar form of the if-matrix: With 
the equations of motion SS/Sa^f^^ = —{2-k)~^ e^^^^ d^bx + 
{2Ti)-'^e^"'^dya~i'' = 0, then, the Lagrangian 



(i) 



— 1 —1 

27r 2n ^ 

i=l 

1 1 

H eada H ea^^^da 

A-K Air ^-^ 



(0 



can be reduced to 



(5.10) 



£ — > — ebda H tada 

2tt An 

+ ZRebda^"^ + Aea(p)aa(p). 
2tt An 



(5.11) 



This is identical to the effective action obtained from the 
flux attachment composite particle approach, where &^ 
plays the role of the vortex gauge field whereas is the 
statistical CS gauge field that transmutes the statistics 
of the electrons. 



B. Parton Construction for Topological Insulators 

in D = 3 + 1 



We can formally repeat the parton construction in 
Z? = 3 -t- 1. As before, we postulate that electrons are 
fractionalized, are consists of p -1-1 partons, and each par- 
ton is in its topological insulator phase. For each parton, 
we can apply functional bosonization to derive its hydro- 
dynamic theory. Solving the constraints among parton 
densities, we will arrive at multi-component BF theories 
with a topological term (axion or 0-term in 13 = 3-1-1). In 
this section, we assume partons are in a D = 3 + 1 dimen- 
sional topological insulator phase in symmetry class AIII 
or Dili characterized by an integer topological invariant. 
Sec Rcfs. |9l|-p^ for previous studies of time-reversal sym- 
metric fractional topological insulators in Z3 = 3 + 1 in 
terms of the parton construction. 



We thus write down the following L agra ngian C 
ELi ^'•'^ + ^'-''^^^ for partons [see Eq. (|3J0| )], 



+ ^e^-'^'d^a^J^dxa^;^ ~ e.J«A^ (5.12) 
for z — 1, . . . p-th flavors, and 

m 

^ m m 

+ E 9,aid,a'^ ep+i 5] (5.13) 



1=1 



7=1 



for the p -I- 1-th flavor. Here, e = J2i = 1, the parton 
densities are written in terms of the two-form gauge fields 
6^*2 and 6^^ as j^'^ = e'^''^Pd^b{p, etc., and are subject 
to the constraint 

m 

^p.Ap^^^W ^ ,A.-Apa^5(P+i) = J2 e^'^'dXp, (5.14) 

1=1 

(i = 1, . . . , p). Sol ving the constraint, following the simi- 
lar steps as in Eq. (5.6), the resulting effective field theory 
is 



L=-J2e'^^'X^dx{a' + Y: 
I 



87r2 



(5.15) 



For a simple case where m = 1 and p = k — 1, in- 
troducing b^^ := bjj^^ and labeling the k gauge fields 

as a^"^'^'"'''^ = {a\}\a'"^\ . . . ,ajj"^), the effective La- 
grangian is given by 

k 

C^-e'^-^Pb^.dxY.a; 



8tt' 



OX 



(5.16) 



As before in the case of L* = 2 + 1, we can eliminate 
one by one. [The following steps should be compared 
with Eqs. (p^-|5lT|) in the D = 2 + 1 dimensional casc.l 
From 5S/6al = 6S/6b^,^ = 0, 



" -AT] =0 



(5.17) 
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If we assume = const, and neglect monopole configura- 
tions, the first equation gives e^'^^'^d^hxp = 0. However, 
if we take into account monopoles, we claim, from the 
first equation 

bxp^^edxa^, a = l,...,k. (5.18) 

Then ^ = ■ ■ ■ = = a, and hence the La- 
grangian, after eliminating aj^, is 



OTT 



(5.19) 

Further integrating over by using SS/Sb^^ = 0, 

e^^-'^-dx (ka, - AD = 0, (5.20) 
we then arrive at 



r — Mt^ApO 4OXQ ACX 



(5.21) 



[This step should be compared with Eq. ( 3.15 ). 



C. Parton Construction for the Quantum Spin Hall 
Insulator in = 2 + 1 



There are two natural ways to apply the parton con- 
struction for D = 2 + 1 time-reversal invariant insula- 
tors: (i) create the necessary copies of the gauge fields 
in D = 4 + 1 and then perform dimensional reduction 
twice (ii) perform dimensional reduction from a non- 
fractionalized theory in D = 4+1 and then create replicas 
of the relevant fields. We will take the former approach. 
Our result, in f act, matches what would be found by 
taking Eq. (5.21) and simply performing dimensional re- 



duction on A';^ 



We begin with Eq. (4^) and add replicas for the 



(f>,ip,an, g^i^jUpi, and to find the Lagrangians 



have 



^fiiyX 



1=1 



HvX 



^ — ' ZTT 



7=1 



7=1 



HvX 



7=1 



(5.23) 



1 the con- 



For the simple case when m = 1, p = 
straint for the parton densities identifies all of the v[l^ = 



v',9liy = gp.u and Uul 



tion of motion 
SS 



5a 



fj,i,. We can then use the equa- 
'-' ^ (5.24) 



to eliminate a^, and arrive at 



c = 



^fl!^X 

"2^ 



—e^^^A'ifdMx, 



(5.25) 



where fix = (j)'^^^ dxtp'^''^ which is the same for all i from 
the equation of motion 6S/Sa\]l^ = 0. Using SS/Scf)^'^'^ ~ 
and SS/5tp^^^ = we can show that i/''^') = tp and (Z)*^*^ = cj) 
are identical for all i respectively. Eliminating ip and 4> 
we finally arrive at 



(5.26) 



Physically this term implies that on the edge of a frac- 
tional quantum spin Hall system there will be fractional 
multiples of e/2 charge i.e. Q = e/(2k) on an anti-phase 
magnetic domain wall. 



_!!::^(^w_e,^-)a^„« 

ZTT 

(5.22) 

where i = 1, . . . , p. Note that we have crucially included 
the Ci charges in front of the scalar fields 0*"^, x"^^ as they 
arise from the dimensionally reduced A™ , A^^ fields re- 
spectively and would enter the D = 4 + 1 Lagrangian 
with the corresponding charge. For the p + 1-st flavor we 



VI. DISCUSSION 

An effective Chern-Simons field theory approach is 
one of the most successful theoretical frameworks of the 
_D = 2 -f- 1 fractional quantum Hall effect. In this pa- 
per, we aimed to extend it this type of hydrodynamic 
formulation to a broader class of non-interacting, as well 
as interacting topological insulators, in arbitrary dimen- 
sions with a suitable set of discrete symmetries. 

We close with a couple of comments. First, our ap- 
proach relies crucially on the presence of a U(l) gauge 
symmetry. While we have focused on topological insu- 
lators with the electromagnetic charge U(l) symmetry. 
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there are topological phases that preserve non-Abclian 
symmetry (such as spin-singlet topological superconduc- 
tors in _D = 2 + 1 and 3 + 1 that preserve spin SU(2) ro- 
tation symmetry) . Our approach can easily be extended 
to such situations (See Appendix^. There are also var- 
ious topological phases that do not conserve any quan- 
tum numbers, except possibly energy and momentum. 
In particular, these include topological superconductors 
in symmetry class D (Z? = 2) and Dili {D = 3). For 
such topological superconductors, Ref. ^ proposed a BF 
type topological field theory with fermionic degrees of 
freedom. 

Second, we note that in various situations, it may be 
useful to initially consider more U(l) charges than arc 
actually conserved. For example, in the quantum spin 
Hall effect, while the electromagnetic U(l) symmetry is 
strictly conserved, spin rotation symmetry is not. It is, 
however, still a useful starting point to consider a system 
with U(l) symmetries associated with both charge and 
spin rotations around, for example, the z-axis. The sys- 
tem is then invariant under an expanded U(l) x U(l) 
symmetry. Functional bosonization in this case then 
gives rise to a BF theory with Chern-Simons terms. This 
is in fact the usual approach and gives rise to doubled 
Chern-Simons theories with U(l) x U(l) symmetry (more 
generally, with U(l)^ x U(l)^ symmetry). Ultimately, 
however, the quantum spin Hall effect does not depend 
crucially on the presence of the additional U(l) spin sym- 
metry; its stability is guaranteed by a Z2 topological in- 
variant, which has nothing to do with the presence of the 
Sz spin rotation symmetry. Starting from the BF theory 
with U(l) X U(l) symmetries, the "unwanted" Sz conser- 
vation can be broken by introducing monopole processes 
or possibly the Higgs mechanism. The former was dis- 
cussed for the BF theory describing a BCS superconduc- 
tor in Ref. ^ which has a Z2 topological order. It would 
be interesting to see if one started from a U(l) x U(l) 
theory and implemented either symmetry breaking mech- 
anism if one would generate a hydrodynamic field theory 
that is equivalent with what we constructed in this article 
by only using the U(l) charge symmetry. 

Finally, due to the lack of microscopic realizations of 
fractional topological insulators (in particular in D > 
2-1-1), the usefulness of our effective field theory ap- 
proach is not yet entirely clear. It is, however, interest- 
ing to speculate on the possible mechanism of such states. 
In the functional bosonization scheme, we need to "raise 
the level" of the BF term, as it is this term that con- 
trols the ground state degeneracy. In this paper, we have 
explored the parton construction. Another possible way 
to raising the level would be to use a Higgs field. This 
mechanism was discussed in the fractional quantum Hall 
effect, in particular, in non-Abelian fractional quantum 
Hall states. Together with functional bosonization, this 
method may allow us to discuss the same kinds of frac- 
tional states obtained by the .parton construction, as well 
as different kinds of states.L30 
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Appendix A: Non-Abelian functional bosonization 

Here, we describe the functional bosonization for non- 
Abelian currents. We start from the generating func- 
tional ^[^J;''] for the correlation functions of non-Abelian 
currents, where A*^^ is a Lie-algebra valued external field. 
We rewrite ^[^Jf] as 

= f V[a,]l[l[5[a, - A';nz[a,]. (Al) 
The delta functional can be written as 

X fl 

^AfpHJ] n e,,^...^p\{5[r^p[a]-r^p[A^-]l 

X f_L,iy.X,... a 

(A2) 

where fap[a] = Saa^ — dj^Qa + [fla, ap], and Yia runs over 
the Lie algebra index, and ^ runs over n = D[D — 
l)/2 independent directions. The Jacobian AFp[a] can 
be expressed in terms of a functional integral as follows. 
By integrating over A"^^ in Eq. (A2), 



V[A^]\{\{5[a^-A^] 
AfpM j V[A^] 

xH n w-n'^[/"/^w-/"^[^]]- (A3) 



(Here we denote A°^ = to lighten notations). Con- 
sider 

(*) := / V[A]\{f{5{rc.p[a] - /„",[A]) . (A4) 

We expand /a^[^] as 

Up[A] = fo^M + D^[a]5Ap - Dp[a]5A^ (A5) 
where A = a + 5 A. Then, 

a</3 

V[5A^] W W 5{Dl[a]5Ap ~ D^[a]M„). 

X a.fS 

(A6) 
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This can be rewritten with a bosonic auxiUary field j3ap 
as 



(*)= / V[p^p,5A^] 

exp* j rf^x p% {Dl[a]5Ap ~ Dl[a]5A^) , 



a<[j 



(A7) 



where (3ap = —Ppa and Pap = for (3 = a. Tlius 



AFp[a]"i = y" P[/3a/3,M^] expi j d"" x PlpDl[a]5Ap. 

(A8) 

The functional determinant can then be written, in terms 
of fermionic ghosts Cq/j and , 

AFp[a] = J 'D[cai3,Cfj,]cxpi J d^xc'^pDa[a]c'p. (A9) 
Then, the generating functional can be written as 
ZlA"""] = J V[a, b, c, c] exp iS, 
S = Jd^xtr c^^...e^"'-'^'^D^[a]cp 

+ 6^,...e^''-"^ {f^p[a] - + L[a]] , (AlO) 

where we introduced L[a] by Z[a] = expi j d^ xir L[a]. 
Let us now introduce auxiliary fields 

j d^xtr c^^...e^"'■■■°'^D^[a]cp 
+ L[a- h] 

- iilh^h^ - 2x/i^c^)] . (All) 

Here, I and x a-re a bosonic and fermionic scalar, respec- 
tively, and /ip is a bosonic. Lie-algebra valued field, /i^ = 
h"" ^ta- The integration over Z sets h^h^ = /i^ = 
and then we go back to the original action. The action 
is invariant under the following BRST transformation: 

6a^, = c^, 6h^ = Cp, 

5x = 1, 61^ 0, 

5cf_, = 0, (5c^i.... = -26^1...., 

= 0. (A12) 

It is possible to write the action as 

iS = SG, 

where G = — 



d^xtr [c 



f^Upla]]. (A13) 



Appendix B: Quantization of the BF theory with 
the Chern-Simons and axion term 



From the bosonized action, we can read off, e.g., the 
current-current correlation functions. If we take the BF- 



CS theory (3.4), however, we do not have charge fluc- 
tuations in the bulk of the system; with only terms of 
topological origin (such as BF and CS terms), there is no 
dynamics. This situation corresponds to the limit where 
the band gap (mass) is taken to be infinitely large. To see 
charge fluctuations, we need either to have other terms 
than BF and CS terms, or to introduce a boundary to a 
system. The former means we keep the band gap finite, 
or consider electron-electron interactions. Let us thus 
consider the Lagrangian with the Maxwell term 



JC = -^6^e^'^^9.(aA - ^a) + ^e^^^^^O^a^ 
47r 47r 



1 



Airg 



2 - M 



(Bl) 



where is a coupling constant, A^ is the dynamical elec- 
tromagnetic U(l) gauge field and F^i, = d^A^, — d^Afj^ is 
the field strength. [A^ is not a static source as in Eq. 
( |3.4|) ]. Given the equation of motion SS/5bfj = 
Ffj,^[A] = ff,4a], the MaxweU term -{i-Kg^)-^ Ff.^F^"' 
gives rise to the Maxwell term for a^, — (47rg^)~^/^,y/^''. 
Alternatively, integrating over A^^ gives rise to a current- 
current interaction cx j^{x)K^i,{x — y)j'^{y) where ac- 
cording to the bosonization rule, 7'^ cx e^^'^^d^bx. 

The BF-CS-Maxwell theory (pl^ can be canonically 
quantized by going through the Dirac quantization pro- 
cedure. Some non-zero canonical commutation relations 
are given by 



Hx)M^')\ 



~i27rCh 



[a,(f),^^(f')] = [A,(a-),^i(f')] =i5l5^^\x-x'), 

(B2) 

where x = (x^, x"^) is the spatial coordinates, i, j run over 
the spatial components of the coordinates i,j = x,y, and 
TT* and it\ are the canonical momentum for a; and Ai and 
are given by 



(B3) 



where Ei = d{)Ai — diA^. The canonical commutation 
relations between current operators are given by 

[/(x),/(a-')] =0, 

[j\x),f{x')]=^^d,5^^\x--x-'), 





An ■' 


^4^^"«- 




+ ^e'^bj - 


ng^ 



An 



(B4) 
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Similarly, in D = (3 + 1) dimensions, one can use the 
bosonized action ( 3.10 ) to read off the current-current 
commutation relations in the bulk, once we allow charges 
to fluctuate, by adding an Maxwell term, say. Let us 
consider the Lagrangian, 



+ ^^''"^'d.^a^dxap - -l^F^^F'^'', (B5) 



where as before is a dynamical U(l) gauge field. 



The BF-axion-Maxwell theory ( p35| ) can be canonically 
quantized by going through the Dirac quantization proce- 
dure. As one may infer from the (half-integral) quantum 
Hall effect at the surface of = 3 -I- 1 dimensional topo- 
logical insulators, the current-current commutators are 
"anomalous" only at the location where 9{x) changes, 
i.e., at an interface between topologically trivial and 
topologically non-trivial insulators. For simplicity, we 
assume 9(x) depends only on the (z) component of 
the spatial coordinates. Some non-zero canonical com- 



mutation relations are given by 

[6™3(f),fo„3(f')] = ^^6""<5(3)(.f-.f' ), 

[a™(f),^^(f')] = [Am{x),7rl{x')] ^ i5l,5^^\x ~ x'), 

(B6) 

where to, n run over the spatial components of the co- 
ordinates m,n = x,y, and tt* and Tr\ are the canonical 
momentum for a, and At and are given by 



1- ' 



:E,,. 



(B7) 



The canonical commutation relations between current 
operators are given by 



167r3 
647r4 ' 



(B8) 



This is exactly the same as the quantum Hall effect in 
the topological insulator in 2 -I- 1 dimensions, but with 
the integer Ch replaced by 83,6. 



^ The Quantum Hall Effect, edited by R. E. Prange and S. 

M. Girvin (Springer, New York, 1987). 
2 M. Z. Hasan, and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

^ X.-L. Qi, and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011) . 

* C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005); C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 
226801 (2005). 

^ R. Roy, Phys. Rev. B 79, 195321 (2009). 

^ B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 

106802 (2006). 

J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) 
(2007). 

* R. Roy, Phys. Rev. B 79, 195322 (2009). 

R. Roy, Phys. Rev. B 79, 195321 (May 2009) 
^ L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 

106803 (2007). 

^" L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
" X.-L. Qi, T. Hughes, and S.-C. Zhang, Phys. Rev. B 78, 
195424 (2008). 

See, Y. Wada et al, Phys. Rev. B 78, 214516 (2008), S. 
Murakawa et al., Phys. Rev. Lett. 103, 155301 (2009), S. 
Murakawa et al, J. Phys. See. Jpn. 80, 013602 (2011), and 
references therein. 

A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 
wig, Phys. Rev. B 78, 195125 (2008). 
^'^ S. Ryu, A. Schnyder, A. Furusaki and A. W. W. Ludwig, 
New J. Phys. 12, 065010 (2010). 



A. Yu Kitaev, AIP Conf. Proc. 1134, 22 (2009). 

A. M. Essin, J. E. Moore, and D. Vanderbilt, Phys. Rev. 

Lett. 102, 146805 (2009). 

S. Ryu, J. E. Moore, and A. W. W. Ludwig, Phys. Rev. B 
85, 045104 (2012). 
1* Z. Wang, X.-L. Qi, and S.-C. Zhang, Phys. Rev. B 84, 
014527 (2011). 

K. Nomura, S. Ryu, A. Furusaki, and N. Nagaosa, Phys. 
Rev. Lett. 108 (2012) 026802. 
2° Q. Niu, D. J. Thouless and Y.-S. Wu, Phys. Rev. B 31, 
3372 (1985). 

2' X.-G. Wen and Q. Niu, Phys. Rev. B 41, 9377 (1990). 

X.-G. Wen and A. Zee, Phys. Rev. B 46, 2290 (1992). 
2^ X.-G. Wen, Adv. Phys. 44, 405 (1995). 
^'^ G. Y. Cho and J. E. Moore, Annals Phys. 326 1515-1535 

(2011). 

P. Nikolic, arXiv: 1108.5388; P. Nikohc, T. Duric 
and Z. Tesanovic, arXiv: 1109 . 0017; P. Nikolic, 
arXiv: 1206. 1055. 

A. Vishwanath, and T. Senthil, arXiv: 1209 . 3058. 

M. C. Diamantini, and C. A. Trugenberger, 

arXiv:1112.3281. 

S. C. Zhang, T. H. Hansson, and S. Kivelson, Phys. Rev. 
Lett. 62, 82 (1989). 

J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 
^° J. K. Jain, Phys. Rev. B 40, 8079 (1989). 

A. Lopez and E. Fradkin, Phys. Rev. B 44, 5246 (1991). 
E. Fradkin, Field Theories of Condensed Matter Physics, 
2nd edition, Cambridge University Press (Cambridge, UK) 



18 



(2012). 

B. Blok and X.-G. Wen, Phys. Rev. B, 42 8133 (1990). 
34 B. Blok and X.-G. Wen, Phys. Rev. B, 42 8145 (1990). 

3^ D. C. Mattis and E. H. Lieb, J. Math. Phys. 6, 304 (1965). 
3^ A. Luther and I. Peschel, Phys. Rev. B 9, 2911 (1974). 
3^ A. Luther and V. J. Emery, Phys. Rev. Lett. 33, 589 
(1974). 

S. Coleman, Phys. Rev. D 11, 2088 (1975). 
3^ S. Mandelstam, Phys. Rev. D 11, 3026 (1975). 
4° F. D. M. Haldane, J. Phys. C: Solid State Phys. 14, 2585 

(1981) 

R. E. Gamboa Saravf, M. A. Muschietti, F. A. Schapos- 
nik, and J. E. Solomin, Ann. Phys 157, 360 (1984), and 
references therein. 

A. M. Polyakov and P. B. Wiegmann, Phys. Lett. B 131, 
121 (1983); ibid 141, 223 (1984). 
''^ Bosonization, edited by M. Stone, World Scientific (Singa- 
pore), (1994). 

A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems, Cambridge 
University Press (Cambridge, UK) (1998). 

E. Witten, Comm. Math. Phys. 92, 455 (1984). 

■"^ K. Fujikawa, Phys. Rev. Lett. 42, 1195 (1979); ibid Phys. 
Rev. D 21, 2848 (1980). 

C. M. Naon, Phys. Rev. D 31, 2035 (1985). 

There have been many attempts at formulating a fermion- 
bose mapping for massless relativistic fermions in Z) > 2 
space-time dimensions. A notabte example is the "tomo- 
graphic" bosonization of Lutheic3 whose bosonic action is 
non-local. However, a variant of this approach has been 
used to bosonize densa-^ermi systems as a theory of a fluc- 
tuating Fermi surface.tj 
A. Luther, Phys. Rev. B 19, 320 (1979). 

F. D. M. Haldane, Proceedings of the International School 
of Physics "Enrico Fermi," course 121, Varenna, 1992, J. 
R. Sciirieffer and R. Broglia editors, North-Holland (NY) 

(1994) ; A. Houghton and J. B. Marston, Phys. Rev. B 48, 
7790 (1993); A. Houghton, H. J. Kwon and J. B. Marston, 
Adv. Phys. 49, 141 (2000); A. H. Castro Neto and E. Frad- 
kin, Phys. Rev. Lett 72, 1393 (1994); ibidPhys. Rev. B 49, 
10877 (1994); ibtd51, 4084 (1995). 

E. Fradkin and F. A. Schaposnik, Phys. Lett. B 338, 253 

(1995) . 

^'^ F. A. Schaposnik, Phys. Lett. B 356, 39 (1995). 

^3 D. G. Barci, C. A. Linhares, A. F. De Queiroz, and ,L F. 

Medeiros Neto, Int. J. Mod. Phys. A 15, 4655 (2000). 
^'^ D. G. Barci and L. E. Oxman, Nucl. Phys. B 580, 721 

(2000). 

K. Shizuya, Phys. Rev. B 63, 245301 (2001). 
'^^ F. A. Schaposnik, AIP Conf. Proc. 419, 151 (1998); 
hep-th/9705186. 

A. P. Schnyder, S. Ryu, and A. W. W. Ludwig, Phys. Rev. 
Lett. 102, 196804 (2009). 

A. Ahland and M. R. Zirnbauer, Phys. Rev. B 55 1142 
(1997). 

R. Jackiw and J. R. Schrieffer, Nucl. Phys. B 190, 253 
(1981); J. Goldstone and F. Wilczek, Phys. Rev. Lett. 47, 
986 (1981). 

'5° L. Fidkowski and A. Kitaev, Phys. Rev. B 81, 134509 

(2010) . 

L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 

(2011) . 

^'^ A. Lopez and E. Fradkin, Phys. Rev. B 59, 15323 (1999). 
P. Hosur, S. Ryu, and A. Vishwanath, Phys. Rev. B, 81, 



045120 (2010). 

A. P. Balachandran, and P. Teotonio-Sobrinho, Int. J. 
Mod. Phys. A 8, 723-752 (1993). 

®^ M. Bergeron, G. W. Semenoff, and R. J. Szabo, Nucl. Phys. 
B 437, 695-722 (1995). 

T. H. Hansson, Vadim Oganesyan, and S. L. Sondhi, An- 
nals of Physics 313, 497 (2004). 
^"^ D.-H. Lee and S. A. Kivelson, Phys. Rev. B 39, 024506 
(2003). 

R. Shindou, K. Imura, and M. Ogata, Phys. Rev. B 74, 
245107 (2006). 

M. E. Peskin, Annals of Physics 113, 122 (1978). 
™ C. Dasgupta and B. I. Halperin, Phys. Rev. Lett. 47, 1556 
(1981). 

M. P. A. Fisher and D. -H. Lee, Phys. Rev. B 39, 2756 

(1989); D. -H. Lee and M. P. A. Fisher, Phys. Rev. Lett. 

63, 903 (1989). 
^2 F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 
^3 R. J. Szabo, Nucl. Phys. B 531 525,(1998). 

R. J. Szabo, Ann. Phys. 280, 163 (2000). 

E. Witten, Phys. Lett. B 86, 283 (1979). 

X. -L. Qi, R. Li, J. Zang, S. -C. Zhang, Science 323, 1184 

(2009). 

H. Kleinert, Multivalued Fields m Condensed Matter, Elec- 
tromagnetism, and Gravitation, World Scientific (Singa- 
pore), (2008). 

^® F. Quevedo, and C. A. Trugenberger, Nucl. Phys. B, 501, 
143 (1997). 

B. A. Bernevig, and S. -C. Zhang, Phys. Rev. Lett. 96, 
106802 (2006). 

M. Levin, and A. Stern, Phys. Rev. Lett. 103, 196803 
(2009). 

M. Levin, and A. Stern, Phys. Rev. B 86, 115131 (2012). 
T. Neupert, L. Santos, S. Ryu, C. Chamon, and C. Mudry, 
Phys. Rev. B 84, 165107 (2011). 
*3 L. Santos, T. Neupert, S. Ryu, C. Chamon, and C. Mudry, 
Phys. Rev. B 84, 165138 (2011). 

Y. -M. Lu, and A. Vishwanath, Phys. Rev. B 86, 125119 
(2012). 

®^ D.-H. Lee and S.-C. Zhang, Phys. Rev. Lett. 66, 1220 
(1991). 

S. -C. Zhang, Int. J. Mod. Phys. B 6, 25 (1992). 
A. Zee, Quantum Hall Fluids, in Field Theory, Topology 
and Condensed Matter Physics, edited by H. B. Geyer 
(Springer, Berlin, 1995). 
®* X. -G. Wen, Int. J. Mod. Phys. B 4, 239 (1990). 
X. -G. Wen, Int. J. Mod. Phys. B 6, 1711 (1992). 
X. -G. Wen, Quantum Field Theory of Many-Body Sys- 
tems, Oxford University Press, (2004). 
J. Maciejko, X. -L. Qi, A. Karch, and S. -C. Zhang, Phys. 
Rev. Lett. 105, 246809 (2010). 

J. Maciejko, X. -L. Qi, A. Karch, and S. -C. Zhang, Mod- 
els of three-dimensional fractional topological insulators, 
arXiv:1111.6816. 
^3 B. Swingle, M. Barkeshh, J. McGreevy, and T. Senthil, 
Phys. Rev. B 83, 195139 (2011). 

T. H. Hansson, A. Karlhede, M. Sato, arXiv : 1105 . 5031, 
E. Fradkin, C. Nayak, K. Schoutens, Nucl. Phys. B 546, 
711 (1999). 

E. Fradkin, C. Nayak, A. Tsvelik, F. Wilczek, Nucl. Phys. 
B 516 704, (1998). 



